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Cherenkov images are multi-fractal in nature. We show that multi-fractal behaviour of Cherenkov images 
arises due to multiplicative nature of pair production and bremsstrahlung processes in the longitudinal shower 
development passage. 

1. Introduction: 

Extensive Air Showers (EAS) [1] are produced when VHE/UHE primary photons ,protons and other high 
Z nuclei enter atmosphere from the top and produce Cherenkov radiation which can be studied by imaging 
technique or by measuring lateral distribution. Cherenkov images recorded by a TACTIC-like 7-ray telescope 
represents the two dimensional distribution of Cherenkov light pattern produced in the atmnosphere. In the 
Hillas parameterization approach, image is approximated as a ellipse and parameters like image "shape" and 
"orientation" are calcualted. Hillas parameters like Length, Width Miss, Alpha, Distance ,Size etc are are 
essentially a set of second order moments. Among these parameters 'a ' and 'Width' are known to be pow- 
erful parameters for segregation of 7- ACE from cosmic-ACE background. Using this approach, present day 
telescopes have been able to reject cosmic ray background events at 98 % level, while retaining up to 50% 
of the 7-ray events from a point source. Hillas parameters are found to be good classifiers for small images( 
close to telescope threshold energy) but tend to fail for large images (of higher primary energy of the order of 
10 Tev or higher) as too many tail pixels are part of the image. Again it is diffcult to study fainter compact 
sources and 7-rays coming non-compact sources or of a diffusive origin. Up to 50% loss of actual image due 
to parameterization is also a big constraint. 

It is interesting to observe that various products of EAS, like Cherenkov photons and particle density distribu- 
tions ( on the ground) have shown multifractal features. Multifractal nature of density fluctuations of particles 
has been experimentally verified and Lipshitz-Holder exponent distribution of EAS has been found to be sensi- 
tive parameter to indentify the nature of individual EAS [2]. We have shown that high energy 7-ray simulated 
images of TACTIC [3] like 7-ray telescope are multi-fractal in nature [4,5]. Multifracal nature of Cherenkov 
images have been confirmed experimentally [6,7]. Because of multifractal nature of EAS products, it becomes 
important to understand the basic multifractal process that takes place during the longitudinal passage of EAS. 

Fractals are self-similar objects which look same on many different scales of observations. Fractals are defined 
in terms of Hausdroff-Bescovitch dimensions. Frcatal dimensions characterize the geometric support of a 
structure but can not provide any information about a possible distribution or a probability that may be part of a 
given structure. This problem has been solved by defining an inifinite set of dimensions known as generalised 
dimensions which are achieved by dividing the object under study into pieces ,each piece is labelled by an 
index i=0,l,2....N. If we associate a probability pi with each piece of size U than partition function can be 
defined as 



r(g,r) = limT(q,T,l) 



(1) 



where 
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For unique function r(q) it has been shown that 



T(q,T) = oo 



(3) 



forr < r(q) 



r(g,T) = 



(4) 



for r > r(g). 

This permits to define generalised dimension D q 



(q - l)D q = r(q) 



(5) 



Here q is a parameter which can take all values between -oo to oo. This formalism is called as multi-fractal 
formalism which characterizes both the geometry of a given structure and the probability measure associated 
with it. 

Cherenokov images arise due to EAS and have been found to be multifractal in nature. In this paper we will 
attempt to understand underlying process responsible for multifractal behaviour of Cherenkov images. 

2. Shower development as a multifractal process 

A ultra relativistic 7-ray enters atmosphere from the top and interacts with air molecule to produce electron- 
photon cascade . The radiation length(x) for pair production and bremsstrahlung process is equal in UHE/VHE 
region [1,10]. Particle -photon cascade in the atmosphere is sustained alternately by electron (e~) -positron 
(e + ) pair production and bremsstrahlung process till the average energy per particle reaches critical value E c 
below which energy loss process is mainly dominated by ionization process. 

Shower development can be visualized a process in which a 7-ray of energy E (= 1 Tev (say)) after traveling 
distance x (on average) produces electron-positron pair each having energy Since energy is getting divided 
into two equal parts , we can attribute a resolution of energy E =2 _1 . In the next radiation length both electron 
and positron lose half of their energy (on average) and each radiates one photon. Thus in this radiation length 
there are two photons and two particles (e _1 and e +1 ) each having energy ^ which can be attributed to energy 
resolution E= 2~ 2 . At this stage fraction of photons ( f>i =5 ) is same as fraction of particles ( pi= 5). As 
the shower develops into next radiation length both electron and positron lose half of their energy and produce 
one photon each. At the same time two photons produced in the previous radiation length interact with air 
molecule to produce electron (e~) -positron (e + ) pair. Thus in third radiation length there are two photons 
and six particles each having energy . This stage can be attributed to the energy resolution E=2~ 3 . It 
is important to note that, at this stage the process of equal division in energy between each photon and each 
particle continues but the process of unequal measure between photons and particles begins. In this radiation 
length the fraction of photons is Pi=\ and fraction of particles (electron +positrons) is P2=|- The fourth 
radiation length corresponds to energy resolution E=2~ 4 as total of 16 photons and particles are produced each 
having energy However, there are 10 charged particles (5 positrons + 5 electrons ) and 6 photons, a case 
of unequal energy division and unequal fraction. 

At a distance of nx , the total number of particles and photons is 2™ , each having average energy |£ and on an 
average shower consists of fraction of | particles and | photons even at nth stage. This corresponds to energy 

resolution of E=2~™. At nth stage each particle or photon can be labelled sequentially with i=0,l,2, The 

probability or fraction of particles and photons can be written as Pi=p\p^ k ' ■ The partition function for finite 
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energy can be written as 

r( 5 ,^) = £« (6) 

where Nk is the number of particles and photons. For the simplicity of calculations, we assume that incoming 
energy E is equal to unit energy. This will not make any difference to actual results but integrate it with 
other classical examples of multifractal behaviour, e.g. curdling of cantor set [9]. On the average shower 
development process has a recursive structure similar to cantor set [9] because both processes are inherently 
binomial multiplicative in nature. For T=r(q),we have 

T(q,T(q),E) = l. (7) 

In the limit of E — > 0, the most dominant contribution to this partition function will survive when T=r(q), 
where r(q) is the solution of the equation 

ip\E- T U> +p\E- r ^) n = 1 (8) 

Above equation can be easily solved to get r(q). 



3. Case of equal energy and unequal fractions 

In each radiation length the measure of photons and particles fluctuates but on the average the shower consists 
of | positrons and electrons and | photons. In the present case E= |, pi= |, p 2 =\. Using above equation and 
Stirling approximation, we have 

I ln(p\+pl) 

Dq q-1 InE (9) 

For q=0, D gives fractal dimension. For q=l, Di is the information dimension which encodes the entropy 
scaling and q=2, D 2 is the correlation dimension which measures scaling of two point density correlation. 
Apart from Dq,Di ,D 2 there are infinite set of other exponents from which information can be obtained by 
constructing an equivalent picture of the system in terms of scaling indices 'a' for the probability measure 
defined on a support of fractal dimension f(a). f(a(q)) is the fractal dimension of the set. In case all D q 's are 
equal to fractal dimensions Dq than f(a(q)) collapses to a point f= a = Do indicating that scaling behaviour 
of a fractal measure is of single type. This is situation which involves neither probability variation nor length 
variation. 



4. Case of unequal energy and unequal fractions 

It is possible that as the shower develops deep into atmosphere, division of energy may not remain equal. 
With the possibility of E\ ^ E 2 ^ \ and/or Pi ^= \ and p 2 ^ § shower development process can still be 
described as a multifractal process because either the energy or the probability or both need to be different for 
a multifractal process. The partition function for this case can be written as 

{p\E- T[q) +P q 2 E- T{q) Y = l. (10) 

where E\ + E 2 =E and pi + p 2 =p. Above equation can be numerically solved for r(q) if values of Ei,E 2 ,pi 
and p 2 are known and fixed for each radiation length. 
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Figure 1. Variation in the muon (4-layer) rate over 

5. Case of loss of energy 

As the shower develops deep into atmosphere most realistic situation is that there will be energy loss in all 
radiation lengths. So a generalized scenario can be described in which E\ ^ E 2 ^ \ and p\ ^ i and P2 ^ § 
with the condition E\ + E 2 < E in each radiation length. Each shower is unique because shower development 
process is random and there are no unique values of E\, E 2 , pi and p 2 . 

6. Simulation studies 

Any simulation / experimental arrangement to measure EAS will consist of set of detectors which can measure 
only a sample or fraction of EAS products. This sample may consist of distribution of charged particles / 
photons. The set of detectors which measure this distribution represents a geometric support and multifractal 
measures can be related to this geometric support. This can be done by calculating multifractal moments 
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of this distribution and obtain generalized dimensions from scaling properties of multifractal moments. The 
definitions of multifractal properties given in previous section are not defined with respect to any support and 
hence can not be applied directly to simulated / experimental data. 

Simulated cherenkov images were generated for 7-rays and protons using CORSIKA code for TACTIC config- 
uration [4,5 and references therein]. 7-rays of energy 50 TeV and protons of energy 100 TeV are considered for 
multifractal studies. Each simulated image is divided into M= 2 V where ^=2,4,6,8, is the scale. The multifractal 
moments are 

M ,k 



where kj is the number of photoelectrons in the kth cell and N is the total number of photoelectrons in whole 
image. G q shows a power law behaviour with M, i.e. G q = M T ^ where r(q) is related to generalized 
multifractal dimension by 

r(q) 

We have calculated average values of D q for 1000 images each for 7-rays and protons. We have repeated above 
studies for 30 TeV 7-rays and 60 TeV protons. No energy dependence of D q on q was observed. 



D q = -^L (12) 



7. Results: 

As the shower develops deep into the atmosphere, values of E\,E%,p\ and P2 in various radiation lengths may 
not remain fixed but may fluctuate. It is not possible to solve equation (4) for all fluctuating values in each 
radiation length. Hence we have to take average values. Figure 1 depicts D q dependence on q for various 
average values. Continuous curve has been obtained in figure 1 for the case Ei=E 2 =\, Pi=\ and f>2=§ ■ 
Simulated data corresponding to protons is represented by (+) sign and simulated data corresponding to 7-rays 
is represented by (x) sign. D q values of each cherenkov image corresponding to 7-rays and protons were 
calculated for 1000 images and average was computed. In figure 1, these average values of D q for both protons 
and gamma-rays have been shown. It is clear from figure that theoretically known average values of shower 
progression does not match with average values of simulated results. 

Figure 2 shows average values of D q versus q behaviour for five different cases. Curve 'a' is actually part (of 
continuous curve) of figure 1, shown here for the purpose of comparison. Average values of D q corresponding 
to proton (+ sign shown as curve d) and 7-rays ( x sign shown as curve e) have been replotted. The continuous 
corresponding to 'b' and 'c' have been obtained . Curve 'b' corresponds to £7i=||, _E 2 =f§, Pi = | an< ^ P2=f • 
Curves 'b' correspond to 23 % loss of energy in proton initiated showers in all radiation lengths on the average 
and curve 'c' corresponds to average loss of 18 % energy in all radiation lengths by 7-ray initiated shower. 
Curve 'a' corresponds to case of no energy loss. It is clear that theoretical model of shower development 
considered with average loss of energy in each radiation length is a good approximation of showers. 



8. Discussion: 

In this paper we have theoritically calculated generalized dimension D q of shower development process based 
on simple model for the case of average values of energy energy and probability. We have considered the 
case of no loss of energy and average loss of energy per radiation length in the shower development process. 
We have also simulated Cherenkov images corresponding to 7-rays and proton initiated showers. Average 
values of D q have been calculated for both types of images. Comparison of average theoritical values and 
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Figure 2. Variation in the muon (4-layer) rate over 



average simulated values are given in figures 1 and 2. It is clear that theoritical model of shower development 
considered with average loss of energy in each radiation length in this paper is good approximation for showers. 

It is interesting to observe that initiated EAS progression comes close to the mathematical definition of classic 
cantor set. Curves b and c in figure 2, represent two scale fractal measure of cantor set in which a fraction 
of energy at each stage gets removed. Such processes have been discussed [15] in detail [page 86]. The 
identification of EAS model as a multifractal process is actually experiencing a cantor process or cantor dust 
in actual physical world. 

In EAS process there is division of energies in terms of charged particles and photons as the appearences of 
charged particles, in any radiation length, means production of cherenkov radiation and appearence of photons 
in any radiation length means absence of cherenkov radiation. The cherenkov image formed on telescope is ac- 
tually superimposition of cherekov photons produced at different heights. Multifractal measures are associated 
the distribution of data on a geometric support. This support may be ordinary plane surface, photomultiplier 
tube (PMT) camera or a fractal itself. 



Multifractal nature 



7 



For particle detection, it has been shown by Kempa[13] that multifractal structures of density fluctuations 
of charged particles near EAS can be used to distinguish between 7-ray and hadron content of cosmic ray. 
The method of multifractal moment analysis has been applied to simulated data as observed by KASKADE 
experiment to distinguish EAS lateral distributions due to various types of primaries [12]. The fractal moments 
are part of the sample of observables used for KASKADE data anlysis [13]. 

Another important point that is clear from this studies is that so called a 'toy model' of EAS can match detailed 
simulated CORSIKA results just by introducing a small loss term in each radiation lenght. 

The different values of E\, E 2 , Pi and p 2 fitted in figure 2 are not unique values. Many other combinations 
of energies and probabilites will produce same resutls. It is because EAS process is random in nature and all 
combinations of different values of energies and probabilites are equally and uniquely probable. It is important 
to note that there cannot be any unique values ascribed to energies and probabilites in EAS. 
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